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Abstract: We consider the constraints for a curvaton with mass m ~ 0{1) TeV and 
(^ ' show that they are not consistent with a purely quadratic potential. Even if the curvaton 

f~^ ■ self-interactions were very weak, they must be accounted for as they affect the dynamical 

evolution of the curvature perturbation. We show that the only TeV-mass curvaton in- 
teraction potential that yields the correct perturbation amplitude, decays before the dark 
matter freeze-out, and does not give rise to non-Gaussian perturbations that are in conflict 
k>( \ with the present limits, is given by l^nt = o^ jM'^. The decay width of the curvaton should 

5-H ' be in the range F = 10~^^ — 10^^^ GeV. The model typically predicts large non-linearity 

parameters /nl and (^nl that should be observable by the Planck satellite. We also discuss 
various physical possibilities to obtain the required small curvaton decay rate. 
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1. Introduction 

In the curvaton mechanism Q, the primordial perturbations originate from quantum fluc- 
tuations of a light scalar field a which during inflation gives a negligible contribution to 
the total energy density (for a review, see [Q]). However, the predictions of the curvaton 
model are quite sensitive to the form of the curvaton potential and the dynamics before 
its decay. As discussed in HI, §, g, |6|, 0, g, H, even small deviations from the extensively 
studied quadratic potential (see e.g. [0]) can have significant effects. This might appear 
paradoxical. However, one should bear in mind that the curvature perturbation (^ is a small 
number, which in the AA^-formalism [^] is the difference between the e-folds of two sep- 
arate FRW universes. Their expansion history is determined by the field dynamics, which 
in the presence of even small non-linearities can lead to widely differing outcomes. The 
number of e-folds in a given separate universe is a large number and its evolution is smooth; 
however, the difference, which is a very small number, can be subject to highly irregular 
oscillations, as was shown in |^. The situation is even more pronounced when discussing 
the non-Gaussianities, which are determined by the field derivatives of the e-folds; here the 
non-linearity of the curvaton potential is essential |^. 

When decaying, the curvaton does not have to dominate the energy density of the 
universe. It is enough that its perturbations are the dominant ones. However, in such 
a subdominant curvaton case the decay products of the curvaton and the inflaton should 
thermalize so as to convert the initial curvaton isocurvature perturbation into an adiabatic 
one ^. In particular, the curvaton should decay before cold dark matter particles freeze out; 



^This is a real challenge for model building. Both the inflaton and the curvaton must decay into the 
visible sector degrees of freedom. There are only few models where the curvaton carries the Standard Model 
charges; a curvaton |H, H, [13, lisl based on the flat directions of MSSM (for a review, see ][L4|) is an example; 
in few cases the MSSM curvaton may dominate the energy density while decaying ||l3, hs, |l5{. 



otherwise the CDM perturbation would not be adiabatic. Such a requirement translates 
into a constraint on the decay rate F of the curvaton. As was pointed out already in [^, P, 
the magnitude of T is also essential for obtaining the observed amplitude C ~ 10~^: given 
some initial condition for the curvaton field a, for a fixed inflationary Hubble rate H^,, the 
curvaton has to decay at a specific time in order to produce the correct perturbation. 

In the present paper we will first address the general question of whether a TeV mass 
scale curvaton is possible at all, given the constraints on the curvature perturbation am- 
plitude, and on curvaton decay rate. In section we demonstrate how in the specific 
case of a TeV mass curvaton the free field assumption is not consistent with theoretical 
and observational constraints and that in this case curvaton self-interactions cannot be 
neglected. In section y we scan numerically the parameter space and compute the curva- 
ture perturbation, study different curvaton potentials and conclude that only a curvaton 
with a self-interaction of the form oc a^ is feasible. In section |^ we discuss the origin of 
small curvaton decay width, F ~ 10~^^ — 10~^^ GeV. Although our treatment assumes a 
real curvaton field, in the particle physics context (such as in MSSM) fields are typically 
complex, and we point a possible mechanism for generating a small curvaton decay width 
that makes use of the kinematical blocking due to the non-zero curvaton VEV. Finally, 
section contains a discussion of the results and possible suggestions for future work. 

2. Importance of self-interactions 

For a quadratic curvaton potential, both the field a and the perturbation Sa have the same 
equation of motion [0]. Hence the relative field perturbation stays constant until decay, 
5a /a = da-t/a-t] here and throughout the paper we adopt the notation where * denotes the 
initial value. Thus one can write the perturbation as 

C ~ — reff ^ 10-^ , 

(Tit 

where H^/a.^ gives the initial perturbation amplitude in the curvaton, and rgfr is the effi- 
ciency factor that can be approximated quite well by the energy fraction at the curvaton 
decay @]: 

Pa 



TeS ~ ?'dcc = 



Pr + Pa 



decay 



Relating a^ and r^ from ^m?a1/'iM'^^H'^ ~ r^, and noting that rjec < 1, we find the 
constraint on the initial curvaton energy fraction 



1 m? , , 

n<-^TTjT- (2-1) 



ec^Ml 



In the free curvaton case r^ec also determines non-Gaussianity through the simple relation 
UmI /nl = 5/4rdec- Very roughly, observationally |/nl| < 100, which implies the constraint 



The limits (|2.1| ) and (2^) are well known. However, there is more. Since the observed 
perturbations are adiabatic to great accuracy, the curvaton must decay before dark matter 
decouples. For each set of the initial conditions, {H^,,r^), there is a relation between rdec 
and the effective decay constant F given by 



'"dec 



Pa 



Pv + Pa 



(2.3) 



H=r 



Here we assume implicitly a perturbative curvaton decay, characterized by decay width F, 
but F could stand for any effective inverse decay time and thus the following discussion 
should hold, at least roughly, also for a non-perturbative curvaton decay as discussed in 
p!^ (note however that non-perturbative curvaton decay could turn out to be a source of 
a considerable non-Gaussianity ]T^). 

The exact evolution of the energy densities is difficult to solve analytically. However, 
we can approximate the curvaton evolution by dividing it up to three phases: 



1. When V" 



m^ < 



H^, the curvaton is effectively massless, so the field value stays 



constant, a = a^, 



2. When V" = m? > H^ the curvaton oscillates in the quadratic potential, and thus its 



-3 



energy density approximately scales as p^^ oc a 
3. The curvaton oscillates until H = T, whence it decays. 
Solving the Friedmann equation for the regime where m? > H^ then yields 
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Using the above result we can solve eq. (|2.3|) to give 
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(2.4) 



To make sure that the curvaton decays before dark matter decouples, F needs to be large 
enough. For, i.e., a LSP originating from MSSM with a decoupling temperature Tlsp ~ 
0(10) GeV, one would obtain as a conservative limit F > 10^^^ GeV, which we will adopt 
as a benchmark value for the rest of the paper. 

In Fig. |l| we have plotted the conditions ( |2.1[ ), ( |2.2D and ( |2.4D for m = 1 TeV; the first 
two appear as two horisontal dashed lines, and r* is limited to be between those two. The 
condition ( |2.4| ) appears as black solid line and the allowed region is to the right. Thus, 
naively, there would appear to be a large parameter region with r^: ~ 10^^^ . . . 10^^^ and 
H* ^ 10^ GeV where a TeV-mass curvaton is allowed. 

However, since the curvaton must decay, it has to have interactions; integrating these 
out would result in an effective theory with some self- interactions. Hence the curvaton 
potential cannot be purely quadratic, although the self-interaction can be very weak, either 
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Figure 1: Parameter space of the quadratic curvaton. r, must be above the red dashed line to 
produce C, ^ 10~^ (equation (2.1)) and below the blue dashed line to produce small enough /nl 
(equation (^^)). Furthermore, F is constrained from above, and thus only the parameter space 
to the right of the black solid line is allowed (equation (p^)). The green dotted line illustrates 



the equality of the mass term and a possible self-interaction term in the potential (equation (2.6)) 
for n — 4. For smaller values of n the line moves further to the left. To the right of the dotted 
line the self-interaction dominates, and thus practically in all of the allowed parameter space the 
self- interaction must be taken into account. 



in the sense that it is a Planck scale suppressed, or alternatively, if the relevant mass scale 
is lower, the corresponding coupling constant is very small. Let us therefore write 

V{a) = -m'a^ + -— . (2.5) 



2 M^i 

In order for the quadratic assumption to be consistent, we need to require that 
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throughout the evolution. Since the energy density of the quadratic field decreases monoton- 
ously, it is sufficient to apply this requirement only for the initial conditions. Solving for 
r* such that the magnitudes of the quadratic and non-quadratic terms are equal, we find 
the condition 
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3M^,H^ 



m^M, 
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(2.6) 



We have plotted this condition for n = 4 in figure ffl as the green dotted line. To the 
right of it, the non-quadratic term dominates initially. As can be seen in figure ffl, there is 
practically no allowed region in the parameter space where the quadratic assumption would 
even approximately apply. For smaller values of n, the self- interaction becomes important 
even for much smaller values of i?* and r*, and thus, there is no quadratic regime left in 
the parameter space. 



We thus conclude that even if the curvaton self-interactions are very weak, a purely 
quadratic potential would not be a consistent approximation for a mass m ~ 1 TeV; instead, 
the effects of the self-interactions need to be taken into account. As we will next discuss, 
the self-interactions change the dynamics of the curvaton in a significant way. 

3. The self-interacting TeV mass curvaton 

Let us now study the dynamical evolution of the perturbations in the curvaton potential 
( |2.5| ) in more detail. We calculate the amplitude of perturbations and the non-Gaussianity 



parameters using A A^- formalism [ 11 1 . This equals to solving the set of equations 



n+3 

= a + {3H + T)& + m'a + {n + A)—— (3.1) 






p, = -AHp, + Ta^ (3.2) 

SiT^Mli = p, + ^&^ + Via) , (3.3) 

in two different, causally disconnected patches of the universe. Here we have assumed that 
the curvaton decay can be at least qualitatively described by an effective decay constant F 
and that ultimately the curvaton decays to radiation p,-. 

The initial conditions are given by pr = 3Mp^H^ and a = 0. However, the initial value 
of the curvaton is not equal in the two patches. Instead, we assume the values ai = a^: 
and (T2 = (T* + H^:/2tt, where V{a^) = r^p^. The curvature perturbation is then given by 
the difference in e-folds, C = ^^ = -/V(c* + i/*/27r) — N^a^,), where N is the number of 
e- folds evaluated at a fixed H after the curvaton has decayed. We then scan through the 
space of all possible initial conditions (ff* , r=K ) , and at each point adjust F so that the right 
amplitude of the perturbations is achieved. We then calculate /nl and (7nl for each point, 
given by the expressions 

5 N" 25 A^"' 

JNL = 7:ir7pT and ^nl 



6^/2 — ^^^^ 54^/3 ' 

where the prime denotes a derivative with respect to the initial condition. This results in 
an allowed region of the parameter space that is limited by the following factors: 

1. The initial perturbation in the curvaton field is too small, i.e., the ratio H^,/a^, is so 
small, that even if the curvaton is completely dominant, the final perturbations are 
nevertheless too small. 

2. The isocurvature limit: To make sure that the curvaton decays only to adiabatic 
modes, we must require that it decays before the dark matter particles decouple 
from radiation background. As before, we take this to translate into the requirement 
F > 10-17 GeV. 

3. The non-Gaussianity parameters /nl and ^nl rnust be within the observed limits. 



Here we use the same limits as in [P, which are given by — 9 < /nl < HI [IS] and 
-3.5 X 10^ < 5NL < 8.2 X 10^ pO|. 



-9 < f_NL ■= 111, -350000 ■= g_NL < 820000 
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Figure 2: The allowed region in the parameter space for n = 4. Since the values of F and rdoc 
oscillate in the non-quadratic regime, so do /nl and (?nl- Thus the requirement for the non- 
Gaussianity parameters creates the distinctive striping. Note that the units of H^ are Mpi. 



The first condition limits the allowed parameter space mainly from upper left, the isocur- 
vature limit from lower left, and the non-Gaussianity limit from lower right. 

3.1 Different powers of the self-interaction 

The value of n is a priori unconstrained. However, from Q we know, that if n > 6, the 
field evolves smoothly also when in the non-quadratic part of the potential, and begins 
oscillations only when it enters the quadratic regime. Due to the lack of oscillations in the 
non-quadratic regime, there are no enhancement of the perturbation or the non-Gaussianity 
parameters. Indeed, we have scanned through the parameter space for n > 6 numerically 
and found no allowed parameter space left. 

A relatively similar case is n = 0. Here the field oscillates in the non-quadratic regime, 
but these oscillations are so fast that the behaviour in parameter space is still smooth. 
We have scanned the parameter space numerically to confirm that for n = 0, there is no 
allowed parameter space left. 

Thus we are left with 2 < n < 6. As was demonstrated in [^ there are strong 
oscillations in the parameter space, and thus one can expect enhancement of C in the 
non-quadratic regime for these choices of n. However, after performing a systematic scan 
through the parameter space, we are left with the surprising conclusion that n = 2 is 
disallowed for a 1 TeV mass. For n = 4 there is however still some parameter space left. 
These are illustrated as the black points in figure y. 

We thus come to the following conclusion: A curvaton with a 1 TeV mass, and with a 
potential of the form 
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Figure 3: /nl and ^nl plotted against r* for all points producing ^ '^ 10 ^ and for a fixed 
ff, == 10^2 GeV. 



can produce the observed primordial perturbations only if n = 4. 

3.2 Predictions for n = 4 

The black pixels in figure |2| correspond to those initial conditions {H^ ,r^), which produce 
(^ ~ 10~^, while still decaying early enough to avoid isocurvature and producing non- 
Gaussianity which is within the observational limits. The values of /nl and ^nl change 
widely from pixel to pixel, due to the oscillations present in the non-quadratic regime. In 
particular, they will change sign; hence there will always be regions in the parameter space 
where they are sufficiently small. However, the average absolute value of both /nl and 
5fNL is so large, that the values of /nl and 5nl in the allowed regions are, except for some 
singular points, typically detectable with future experiments, i.e., the Planck satellite. This 
is demonstrated in figure ^ for a fixed value of H^: = 10^^ GeV. 

In figure § the corresponding values for F are given, which demonstrates that the al- 
lowed range is given roughly by 10~^^ . . . 10^^^ GeV, i.e., the allowed points mostly populate 
the lower limit for F. 

In general, the numerical results can be summarised to be the following: 



A curvaton with a mass of 1 TeV is limited to the region in the parameter space 
where self-interactions suppressed by the Planck mass become important. 
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Figure 4: The values of F corresponding to the ( ^ 10~^ for n — A. The units of iJ* and F are in 
Mpi. Note that in comparison to figure 0, points with too large /nl or gNL have not been removed 
from this plot. 

• Of all the possible powers of the self- interactions, only the case n = 4 has non- 
negligible allowed parameter space. This corresponds to a self-interaction of the type 

Vint OC 0-*. 

• The typical values of both /nl and (^nl are so large in the allowed patches of the 
parameter space, that they are observable with future CMB observations. 

Even though we have here solved the equations numerically for mass of 1 TeV, we can 
deduce the results also for other masses. Prom the analysis of |P] and 0] it is clear that 
enlarging the mass will also enlarge the allowed regions of the parameters space. Thus for 
a 10 TeV mass, the plot of the allowed points in the parameter space would look otherwise 
similar to figure |2|, but the black stripes would be somewhat thicker. Similarly, for smaller 
masses, i.e. 100 GeV, the allowed region will become smaller. For a 10 GeV mass practically 
no parameter space would be left. 



4. Source of small F: Kinematical blocking 

As demonstrated by figure ^, the allowed values for the decay constant F populate the 
lower bound, i.e., are in the range 10~^^ . . . 10~^ GeV. Although we have not constrained 
ourselves to any single particle physics model, for most realistic scenarios this coupling is 
very small. To alleviate this problem, we suggest possible mechanism to implement such a 
small decay constant. 

After inflation the majority of the energy density of the universe is in the relativistic 
degrees of freedom. We strictly assume that they belong to the visible sector degrees of 
freedom which would thermalize with that of the Standard Model quarks and leptons. 



Those degrees of freedom which are coupled to the curvaton wih obtain VEV dependent 
masses, which for a quartic couphng T^nt ~ h?'a'^(jP' is meff ~ h{a). The decay of the 
curvaton is thus delayed until meff < m, ~ C(l) TeV. To estimate when the curvaton 
decays, one has to follow the evolution of the VEV of the curvaton to find out when 
m ~ /i(cr). 

For a real valued curvaton, the curvaton field oscillates through zero, and thus the 
kinematical blocking is lifted for a short part during each oscillation. In such a case. 



however, instant preheating-type phenomena will typically become important |21], see for 
a review p^. For a complex field the field trajectory cannot pass the origin if there is any 
initial rotational motion [£3|. Strictly speaking, our analysis in the preceding section is 
valid for real fields alone. However, starting from the action for a complex field, 

and inserting an Ansatz ip = \(f)\e yields an effective action of the real modulus field with 
an effective mass 

2 2 n2 

m^g = m — U . 

Hence we may still adopt our analysis in the preceding sections for the real field if the 
rotational motion is not too fast. Equivalently one could write the equations of motion for 
the field (p — )■ (pe (without the Hubble friction term) as 

^2 

(p + m^cp — — = , 

where L is a constant associated with angular momentum. The evolution of the radial 
component of the field is similar to the real field case for large values of (j), but near the 
origin the effective potential rises fast, and the radial part is refiected away from the origin. 
If the time that the field spends near the origin is very small compared to the background 
dynamics (Hubble time) , then the evolution of the complex field will effectively be the same 
as for the real field. This requires small enough L. 

Thus for the decay of a complex field to be kinematically blocked, we require that 
there is sufficient rotational motion as to prevent the field from traversing too close to the 
origin, yet at the same time to be small enough for the analysis of the preceding sections 
to apply. Here we assume that these conditions are met. 

In the quadratic case, the VEV of the curvaton will be redshifted as a~^'^{t) once the 
flat direction commences its oscillations at H{t) ~ m. In a radiation dominated epoch the 
scale factor is a{t) oc H^^''^{t). Thus 

a{H) = a, (-) " . (4.1) 



Relating the initial field value a^ to r^ used in figures @ and Q with the equation r^ = 
^m'^a'^/SMp^H^, we can solve the value of H at the moment the curvaton starts to decay: 



iJdcc = rn 



1 



/EnhMpiH^ 
In contrast, in the non-quadratic case, the VEV starts to decrease when the curvaton 



becomes massive, H < Hgj 



m^ 



2 + SGa^/Mpj, while the field scales as a (x a '^'^, and 



aiH) = a.[—' 



(4.2) 



This continues until the VEV has diminished sufficiently and enters into the quadratic 
regime. The point of transition is defined by g'^^'^oq ~ ^eq/-^pr After this point the VEV 
behaves as in the quadratic case with 



aiH) = a, 



eq 



H 



H, 



eq 



Now again requiring that the curvaton decays when ha ~ m, we get 



m 



(^{Hdec) = 0"cq ( -77^ j 



Solving Hcq from Eq. (4.2), 



(4.3) 



-f^eq — Hs 
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-'^dec — -"eq 



Inserting the expressions for H^q and Ueq, we get 
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, , m2Mjl, \ 6 



Solving 0",, from r* = u^/Mpj /SMp^ii'^, and inserting it into ( |4.5| ), we get 



(4.5) 



a. 



_ m^ \jm^ + 56 r] r| Mp\ 

'^'^'^ ~ 1 5 ? 4 A I 
23 3i2M^i/i3ri2i/,« 



(4.6) 



The curvaton decay takes place when V ~ -H^dec- The perturbative decay rate is de- 
termined by Frf ~ /i^TTi, which is larger than the kinematical blocking, i.e. F^ ^ F ~ 
Uii^c.- Therefore once the kinematical blocking is lifted, the curvaton decays immedi- 



ately H, H, m, |2|. 
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Figure 5: Value of i^doc plotted for H^ = lO^^ GeV and for masses 1000, 100 and 10 GeV. 



The values of i/dec ~ T for some typical values of iJ* and r^ are given in figure ^. The 
time of decay is relatively insensitive to all the other parameters, and strongly decreases as 
the function of the curvaton mass. A relatively small F ~ 10^^^ GeV can be obtained for 
small masses (i.e. m ~ 100 GeV) for a coupling of order unity, h ~ 0{1). The perturbative 
value of Yd is generically quite large for /i ~ 0(1). To acccommodate the small value 
for r ~ 10~^^ — 10~^^ GeV, the curvaton must couple to the SM degrees of freedom very 
weakly, i.e. h ~ 10~^ — lO"^'^. Unfortunately, within MSSM it is not possible to find such 
a curvaton candidate unless the decay of the curvaton is somehow delayed. 

One such attractive possibility may arise if the curvaton instead of decaying directly, 
first fragments into lumps of supersymmetric matter, known as Q-balls pH]. The formation 
of Q-balls within MSSM has been extensively studied analytically p7[ and numerically 



The Q-ball formation is quite robust even if the curvaton oscillates radially [29|, all that 
we require is that the curvaton oscillations on average feel negative pressure. In presence 
of a negative pressure certain sub-Hubble modes of the curvaton perturbations become 
unstable and leads to the formation of Q-balls. These Q-balls do not decay immediately. 



instead they evaporate into pair of fermions from their surface |26], naturally suppressing 
the effective decay rate [^] (for a review, see [|l4|). However the details of the decay rate 
depends on the charge accumulated in the Q-balls ^. 

There are three possibilities emerge within MSSM which can fragment to form Q-balls. 
They are the udd, LLe and QLd monomials [^] . Some of these monomials are lifted by non- 
renormalizable superpotential operators of dimension d = 5, such as HuLLLe, H^Ludd. 
Typically these superpotential operators are hybrid in nature and would never allow a 
stable non-renormalizable 74-term in the potential proportional to the superpotential itself. 
As a result the potential for such flat directions would look exactly like Eq. ( ^.5| ) with 
n = 4. 



^One of the astrophysical signature of fragmentation of any condensate is the production of very high 
amplitude gravity waves which can be detectable by future gravity wave detectors |3C(] . 
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5. Discussion 

In this paper we have considered requisites for the existence of a TeV mass curvaton. Such 
a curvaton has to produce the observed perturbation ampUtude and not too much non- 
Gaussianity; moreover, it has to decay before the CDM freeze-out. These constraints fix 
the range of the initial conditions which turns out to be such that the quadratic term 
in the curvaton potential cannot dominate over possible higher-order terms for the whole 
dynamical range. Hence we conclude that a purely quadratic curvaton potential would not 
be a consistent description of a TeV mass curvaton. However, the presence of non-linearities 
in the field equation of motion then very much modifies the outcome for value of the 
curvature perturbation. Interestingly enough, as discussed in Sect. |3|, we find that the only 
viable curvaton potential that satisfies all the observational constraints is 1^ = rf^a^ jl ■\- 
a^ /M'^. Moreover, the curvaton decay rate should be in the range T = 10~^^ — 10~^^ GeV. 
Note that in the case where the curvaton energy density is subdominant at the time of 
decay, the curvaton does not necessarily have to decay before baryogenesis, which can be 
a process that takes place among the inflaton decay products. However, the decay should 
be able to produce thermal CDM particles so that the CDM perturbation is adiabatic. 

As we point out, MSSM flat directions which are lifted by d = 5 operators would 
produce exactly the form of the potential required for the TeV mass curvaton. Here we 
should emphasize the fact there are directions with no ^-terms so that although the flat 
direction field is complex, its modulus would, to a good approximation, have the same 
equation of motion as the real field considered in Sect. ^. Rotational motion in the complex 
field plane would effectively only modify value of the mass parameter, and as long as the 
angular momentum is small enough, the results of Sect. ^ should apply. However, the 
rotational motion is essential for the kinematical blocking that protects the flat direction 
field from decaying as long as the field is away from the origin. 

As discussed in Sect. Q, kinematical blocking is also modified the presence of non- 
linearites. We find an effective F that is small, but perhaps not small enough. For some 
regions in the parameter space, especially for small curvaton mass parameters, the correct 
perturbation can nevertheless be generated. We should however stress the approximate 
nature of our estimate. A proper calculation would account not only for the actual ro- 
tational motion but also the dynamics (and possible backreaction) of the kinematically 
blocked degrees of freedom. This would provide a major numerical challenge but might 
nevertheless be worth the attempt. Meanwhile, we can only conclude that a TeV mass 
curvaton based on a d = 5 MSSM flat direction remains a possibility. Another natural 
solution may arise if the flat direction fragments to form Q-balls, which then evaporate via 
surface evaporation which can delay the decay rate sufficiently. 

Note also that what really matters is the equation of state, not the time of decay. 
Thus if the curvaton decays too early, the perturbations might still generated if the decay 
products have the equation of state of matter. An example of this could be the MSSM flat 
direction decaying into Q-balls, which would then slowly decay. Another possibility could 
be that the MSSM curvaton would not be reponsible for the amplitude of the perturbations, 
but might add a small component on top of the inflaton perturbation spectra with (very) 
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large values of /nl and ^nl- The production of large non-Gaussianities by the MSSM flat 
directions could in fact be a generic feature that would deserve a closer study. 
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